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The dependence of the ratios of nuclear structure functions is studied by performing 
QCD evolution of nuclear parton distribution functions. The log slope of these ratios is 
very sensitive to the nuclear gluon distribution function. Taking different parametrizations, 
we show that the NMC data on the dependence of i^f rule out the case where nuclear 
shadowing (suppression) of gluons at a; ~ 0.01 is much larger than the shadowing observed 
in the ratio F2 /F2' ■ We also take into account modifications to the DGLAP evolution by 
including gluon fusion terms and see that the effect is small at present energies, and, in any 
case, a strong gluon shadowing is not favored. The region studied {x ~ 0.01) is the most 
relevant for RHIC multiplicities. 



The nuclear parton distribution functions (nPDFs) are an essential ingredient in calculations 
of hard processes involving nuclei. In the framework of the QCD improved parton model it is 
possible to extract the nPDFs from experiments of DIS and to use them for other processes in 
AB collisions. In this framework, a set of nPDFs, or equivalently, nuclear corrections to PDFs 
in nucleons, has been obtaineda'Ba. There are, however, some uncertainties in the determination 
of the initial conditions for the evolution due to the lack of experimental data. For instance, 
gluons and sea quarks are not constrainedu in the EMC region. Also, it is a common belief that 
nuclear gluon distributions are very weakly constrained by experimental data at small x. The 
situation would be similar to that of nucleon PDFs before HERA measured the rapid increase 
of as X ^ 0. However, we have shownll that within the leading twist (LT), lowest order (LO) 
DGLAP framework, the DIS data provide enough constrains to rule out very strongly shadowed 
gluon distributions in the region x>0.01. The key experimental information comes from the 
dependence of the ratios F2'^ /F2 measured by the NMC Collaborations. 



In the case that the dominant contribution to total multiplicities comes from jets or minijets, 
as has been proposed, the measurements of Nch done at RHIC give direct information about 
the initial gluon distributions in nuclei. This has been used in Ref. 13 to parametrize the gluon 
distribution function in the HIJING model. As a result, a strong gluon shadowing has been 
proposed in order to compensate the large multiplicity obtained in the original model. 

In previous works SH we have determined a set of nPDFs (named EKS98) by using available 
data of F2{x, Q^) and Drell-Yan cross sections in nuclear scatterings as constraints. The iterative 
procedure is very similar to that in global fits of parton distributions of the free proton. As a 
result the nPDFs become fixed at an initial scale Qg, and the evolution to > Qq is given by 
the DGLAP equations!. The nuclear case is more complicated as a new variable {A) is present. 
Moreover, as the amount of data does not allow to fix the nuclear ratios for each parton species 

Rtix,Q'^) = "^f.l^^'q'^y 'i = 9, Uv, d^, u, . . . , (1) 

in an unambiguous way, some assumptions were imposed in addition to the constraints from 
the momentum sum rule (that, roughly speaking fixed the amount of antishadowing for gluons) 
and baryon number conservation. The assumptions include the saturation of the shadowing at 
small values of x and the existence of an EMC effect for gluons and sea quarks at large x. 

Notice that the DGLAP evolution imposes some constraints to nuclear gluons, just in the 
same way as it does in the global analysis with nucleons. In particular for small values of x, LO 
DGLAP equations give 

51ogQ2 - 27vr (x, Q^) l^^^^rr, Q ) - i?^,(x, Q )|. (2) 

So, the slope of the evolution of the ratio F2 /F2 is related to the difference of the nuclear 
ratio in gluons and F2 (notice, however, the factor 2x in Rg). As Rp^{x,Q'^) is measured 
experimentally, data on the evolution of the ratios Rp^ constrain Rg once DGLAP evolution 
is taken into account. In particular, the experimental data from the NMC CollaborationB show 
that F2'^ /F2 increases as increases at small values of x. For this reason, R^ cannot be much 
smaller that Rp^ in order for the slope to be positive. We have studied different parametrizations 
for the initial ratios, and shown this fact by performing the DGLAP evolution without making 
use of the approximation @). We will see that the result is in agreement with the argument 
above. 

The initial conditions of the four parametrizations studied can be seen in Fig. 1. Notice that 
the HPCb and the new-HIJINGQ parametrizations are independent. They are used here as 
initial conditions for evolution. HKM0 is based on DGLAP evolution, but the NMC data! on 

dependence of -F^'^/Fg" and DY data have not been taken into account. We think that this 
is the main reason for the difference in gluons and sea quarks as compared with EKS98. 

As can be seen from the figure, all the parametrizations give similar results for quarks 
at the values of x where experimental data exist. However, the largest uncertainties are in 
the gluon distribution function, the extreme case being the new parametrization introduced in 
HIJING which suggests a very strong shadowing. We will make use of this parametrization to 
demonstrate that a very strong shadowing for gluons is ruled out based on the experimental data 
on the experimental data on the dependence of F2^ /F2 ■ For that we perform the DGLAP 
evolution of the nPDFs given by these parametrizations for Sn and C and, then, compute the 
ratios. In Fig. 2 we present the result of this analysis. It is clear that the strong gluon shadowing 
case is in complete disagreement with the data as it presents a negative slope, while all the others 
are positive. 

One could argue, however, that the non-linear terms (expected to be present in evolution 
equations at high enough density) could change this conclusion. In order to check this possibility. 




Figure 1: Initial ratios Rfix^Q^) at Qq — 2.25 GeV^ of Pb, for valence quarks (solid), sea quarks (dashed) and 
gluons (dotted) from the four different parametrizations used. 

we repeat the evolution using the same initial mnditions but now taking into account the 
gluon-fusion terms computed by Mueller and Qiu EJ. In this case, the evolution is modified 
as can be seen in Fig. 3. Including these non-linear terms the slopes become always smaller 
(even more negative in the case of strong gluon shadowing). So we can conclude that using the 
presently available experimental data, a strong gluon shadowing is ruled out at x>0.01 in the 
present framework. 
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FiguKe 2: Comparison of the calculated El and measured dependsfice of the ratio -Ff"/-F^-r 
dataB are shown with statistical errors only. The results for EKSOSHn (solid lines) and HKMO 
dashed) are fgpm the corresponding global DGLAP analyses. The dependence of the HPCl 
and HIJINGLl (dotted) cases is obtained from the DGLAP evolution with initial conditions 

Fig. 1. 
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Figure 3: r^he scale dependence of the ratio jF^ calculated using DGLAP evolution with MQ 
corrections llj, and compared with the NMC dataS. Initial conditions for nuclear effects are taken from 
EKS98 (solid hues), HPC (dashed) and HIJING (dotted) parametrizations. 



